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Abstract 

The Chiral Magnetic Effect (CME) is predicted for mid-peripheral Au-Au collisions 
y^s]v7V=200 GeV at RHIC. However many backgrounds can give signals that make the 
. measurement hard to interpret. The STAR experiment has made a measurement that 

I makes it possible to separate out the background from the signal. An event shape 

analysis is the key for controlling the effects of Transverse Momentum Conservation 
CD ■ (TMC) in the Au-Au colhsions. LHC Pb-Pb cohision ./s^=2.76 TeV has made 

Q I event plane measurements from which we calculate the CME and make predictions 

I about what the event shape analysis should give. 



1 Introduction 



Topological configurations should occur in the hot Quantum Chromodynamic (QCD) vac- 
uum of the Quark-Gluon Plasma (QGP) which can be created in heavy ion collisions. These 
^ ■ topological configurations form domains of local strong parity violation (P-odd domains) in 

0\ , the hot QCD matter through the so-called sphaleron transitions. The domains might be 

detected using the Chiral Magnetic Effect (CME)[1] where the strong external (electrody- 
CN I namic) magnetic field at the early stage of a (non-central) collision, through the sphaleron 

CN ' transitions induces a charge separation along the direction of the magnetic field which is 

perpendicular to the reaction plane. Such an out of plane charge separation, however, varies 
its orientation from event to event, either parallel or anti-parallel to the magnetic field 
(sphaleron or antisphaleron). Also the magnetic field can be up or down with respect to 
^ \ the reaction plane depending if the ions pass in a clockwise or ant i- clockwise manner. Any 

^ ■ P-odd observable will vanish and only the variance of such observable may be detected. 

The STAR collaboration [2] has published a measurement of charge particle azimuthal 
correlations consistent with CME expectations. In Ref.[3] and used by STAR the CME can 
be indirectly approached through a two-particle azimuthal correlation given by 

7= (cos(0i + (/)2-2^if,p)) = (cos(0i-^Rp)cos(02-^i?p))-(sm(0i-^/jp)sm(02-^Rp)), 

(1) 

where "^rp, (pi, 02 denote the azimuthal angles of the reaction plane, produced particle 1, 
and produced particle 2. This two particle azimuthal correlation measures the difference 
between the in plane and out of plane projected azimuthal correlation. If we would rotate 
all events such that \1/_rp = 0.0, then 7 would become 



7 = (cos(0i + 02)) = (cos(0i)cos(02)) " (sm(0i)sm(02)) • 



(2) 



The CME predicts that 7 > for opposite sign-pairs and 7 < for same sign-pairs. There 
are other two particle azimuthal correlation effects that can depend on the reaction plane 
driven by elliptic flow even though the underlying correlation may be independent of the 
reaction plane. These backgrounds are summarized in Ref . [4] . It was pointed in Ref . [4] that 
the difference correlation (5) which is independent of the reaction plane gives a constraint 
on the CME and backgrounds. 

5 = (cos(0i - (f)2)) = {cos{(f)i)cos{(j)2)) + (sm((/>i)sm(02)) ■ (3) 

In Ref. [4] Transverse Momentum Conservation (TMC) is derived and demonstrated that 
if there is no other correlation in the data except elliptic flow TMC will give a negative 
{cos{4>i -\- 4>2 — '2^ Rp)) (7) which is V2 smaller than {cos{4>i — ^2)) (^)- 5 is a negative 
number given by TMC and scales as 1/N (N is the number of particles). Depending on 
other correlations the connection between 7 and 5 could be larger than V2. Also in Ref. [4] 
Local Charge Conservation (LCC) is another important background and the details of LCC 
is found in Ref. [5]. The authors of Ref. [5] point out that for same sign pairs LCC should 
give a small negative sign, while for opposite sign pairs LCC should give a large positive 
correlation. Using the same coupling effect to the reaction plane as TMC, one should expect 
(cos(0i + (j)2 — 2^ijp)+_) (7) equal V2 times (cos(0i — 02)+-) (S) for the LCC. 

Ref.[l] has pointed out that P-odd domains on the surface of the flreball omit same charge 
sign particles in the direction of the magnetic held. The particles that escape the surface 
would be of the same sign while the charge particles moving in the opposite direction would 
be of opposite sign. These particles would run into the flreball and be thermalized and loss 
their direction (quenched). This effect makes it hard to know what the CME should give 
for opposite sign pairs. Also when one considers opposite sign pairs the LCC effect is much 
stronger than the TMC effect making 6 a large positive number. When we scale down by 
the expected V2 to obtain 7 we obtain a larger positive number which is bigger than what is 
measured. 6 (2.0) times V2 (-06) equals 7 (.12) where the measured value of 7 is .06. We need 
some additional negative .06 from some other source. Minijet quenching can increase the 
(sm(0i)sm(02)) term because there is more quenching out of plane than in plane. This could 
account for a reduction of .06. I think it is clear that having three important effects plus 
two quenching effects makes the opposite sign very hard to make an analysis. On the other 
hand LCC, minijets and CME quenching for the same sign can be considered unimportant. 

The paper is organized in the following manner: 

Sec. 1 is the introduction to correlations. Sec. 2 is an analysis of the same sign data. 
Sec. 3 LHC analysis and predictions. Sec. 4 presents the summary and discussion. 

2 Same Sign Analysis 

In performing the same charge sign analysis we will assume that only TMC and CME 
amplitudes are present in the Au-Au y^sjvjv=200 GeV data. We will use (cos(0i-|-02— 2^ijp)) 
(7) and (cos(0i — ^2)) (S) of Ref. [2] in our analysis. Figure 1 is the 7 measurement of 



STAR while Figure 2 is the 6 measurement. We can separate the Transverse Momentum 
Conservation (TMC) effect in the Au-Au colhsions using a STAR event shape analysis. For 
20-40% centrality half of the tracks of the TPC are used in order to determine the reaction 
plane. The other tracks are used to calculate {cos{(f)i + (l)2 — 2^jip)) (7) for same and opposite 
sign. Also {cos2{(j) — "^rp)) {V2) is calculated for these tracks. The V2 measures the shape 
of each subevent. For positive V2 we have more tracks in plane, while for negative V2 we 
will have more tracks out of plane. The value for TMC in the correlator 7 should be zero 
when V2 is zero. The value should be negative when V2 is positive and positive when V2 is 
negative. On the other hand the CME should not depend on the event shape (^2) but only 
on the collision geometry which is given by the centrality. Thus the magnetic field and the 
charge separation is a constant offset. 

In Figure 3 the results are plotted from the STAR analysis which used correlators that 
are proportional to -7. Ref.[4] used measurements from Ref. [2] for centrality range 50% 
to 60% because this should give the best measurement signal to errors. The value for 
(cos(0i + (f)2- 2*Rp)) (7) is -.00045 ± .00005 and {cos{(j)i - 02)) (5) is -.00038 ± .00010. V2 
is 0.06 and from Figure 3 assuming that the data for the same charge sign will be proportional 
to the centrality range 50% to 60%, we can determine the ratio of pure CME to 7 at measured 
(cos2(0_*j?p)) = 0.0 and (cos2(0_*rp)) = .06. Figure 3 gives a ratio of 0.60 ± 0.15. 



2.1 Fit Model 

The fit model that we will use assumes that the three experimental values 7, S and the ratio 
for the centrality range 50% to 60% can be fitted by three equations. 

7 = (cos(0i + 02)) = CME + ScaleT * TMC. (4) 

S = (cos(0i - 02)) = TMC - ScaleM * CME. (5) 
. _ CME _ CME 
ratio - - CME + ScaleT * TMC ' ^ ' 

Where ScaleT is the faction of the TMC which is measure as part of 6 as translated to 7 and 
ScaleM is the faction of the CME which is measure as part of 7 as translated to S. From the 
above equations we have three equations and four unknowns. 



2.2 In Plane and Out of Plane 

Equation 2 and equation 3 show how 7 and 6 can be broken up into an in plane term 
(cos(0i)cos(02)) and an out of plane term (sm(0i)sm(02))- If we would make a simple 
assumption that the in plane (cos(0i)cos(02)) is equal to the TMC value and that the out 
of plane (sm(0i)sm(02)) is equal to the CME value, one can fit the equation 4 and equation 
5 with ScaleT and ScaleM equal 1.0. TMC is equal to -.00041 and the CME is equal to 
-.00004, but the ratio is only .09 which is way off. Let us consider the behavior of the CME 
term. Let us look at the pure CME case. Equation 4 becomes 



(cos(0i)cos(02)) - (sm(0i)sm(02)) = CME, 



(7) 
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Figure 1: (Color) (Solid points) (cos(0i + (1)2 — 2^1' /jp)) (7) calculated for Au-Au ^/sN^=200 
GeV events data acceptance cuts of 0.15 < pt <2 GeV/c and jr^l < 1.0. Event generators HI- 
JING used the same input and acceptance cuts with and without an elliptic flow afterburner. 
Also generators URQMD and MEVSIM where used. Blue symbols mark opposite charge cor- 
relations, and red are same charge. For MEVSIM, HIJING, and URQMD points the true 
reaction plane from the generated event was used for "^Rp. Thick, solid lighter-colored lines 
represent possible non-reaction-plane-dependent contribution from many-particle clusters, 
as estimated by HIJING and discussed in Sec. VIIA of Ref.[2]. 
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Figure 2: (Color) (cos(0i — ^2)) (5) calculated for Au-Au ys/v7v=200 GeV events data 
acceptance cuts of 0.15 < pt < 2 GeV/c and \r]\ < 1.0. We show different charge combinations 
and forward full magnetic field (FF) and reverse full magnetic field (RFF) configurations. 
The data points corresponding to different charge and field configurations are slightly shifted 
in the horizontal direction with respect to each other for clarity. The error bars are statistical. 
Also shown are model predictions described in Sec. VII of Ref.[2]. 
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Figure 3: The charge correlation difference (Aud) — (A-lr) and (A_,_A_)ud — A_|_A_)lr which 
is proportional to —7 for the same sign pairs (squares) and opposite sign pairs (dots) under 
the selection of soft particles low-pt {pt < 2 GeV/c) is plotted as a function of azimuthal 
anisotropy shape of the sub-events. 



while equation 5 becomes 

{cos{4)i)cos{4)2)) + {sin{(j)i)sin{^2)) = -ScaleM * CME. 

(cos((^i)cos(02)) and {sin{(l)i)sin{(f)2)) can be written as 

2(cos(0i)cos(02)) = (CME - ScaleM * CME), 

and 



(8) 



(9) 



2(sm(0i)sm(02)) = -(CME + ScaleM * CME). (10) 

We know that CME is left right symmetric therefore {cos{(j)i)cos{(j)2)) term equation is zero 
thus 

0.0 = (CME - ScaleM * CME), (11) 
or ScaleM — 1.0. This is true in general. 

2.3 Final Fit 

We have seen in the last subsection that in general ScaleM = 1.0 thus equation 4 through 6 
becomes 

7 = {cos{(pi + 02)) = CME + ScaleT * TMC. (12) 
5 = (cos(0i - 02)) = TMC - CME. (13) 
CME CME 



ratio 



7 CME + ScaleT * TMC 
We now have three equations and three unknowns. After doing the fit we obtain Table I. 

Table I. Fit parameters and errors. 



(14) 



Table 1 


parameters 


value 


errors 


CME 


-.0002625 


+ .0000638 -.0000845 


TMC 


-.0006460 


+.0001195 -.0001295 


ScaleT 


.2836 


+.1601 -.1296 



Ref.[4] made a prediction that ScaleT should be equal to V2 or 0.06. However there must 
be some other correlation of same sign charge pairs which makes the effect four times bigger. 
Ref.[6] showed that there is a suppression of near side same sign charge pairs created by the 
strong chromo-electric fields inside the glasma fiux tubes generated by the color glass initial 
state of the Au-Au collisions. This suppression greatly increases the negative TMC effect 
between back to back particles of the same charge sign. 



3 LHC Analysis and Predictions 



At the LHC the ALICE experiment has made measurements Uke STAR but on Pb-Pb at 
yj/vjv=2.76 TeV. We will again performing a same charge sign analysis using (cos(0i + 02 — 
2^Rp)) (7) and (cos(0i — ^2)) (S). Figure 4 is the 7 measurement of ALICE while Figure 5 
is the 5 measurement. Since LHC is at such a higher energy hard scattering or jets play a 
larger role. Let us initially assume that only jets and CME are present in the data. 

Again the fit model will use the centrality range 50% to 60% for the LHC. The value for 
(cos((^i + 02 - (7) is -.0004 and (cos(0i - 02)) (6) is .00125. The V2 value for this 

data is .105. We break up the 7 and 6 into an in plane term (cos(0i)cos(02)) and an out of 
plane term (sm(0i)sm(02))- 

(cos (01 ) cos (02)) = .000425, (15) 

and 

(sm(0i)sm(02)) = .000825. (16) 

3.1 Jets and CME 

Using Ref.[4] as apphed to jets, we can write 

7 = (cos(0i + 02)) = V2Jet. (17) 

S = (cos(0i - 02)) = Jet. (18) 
Thus solving for (cos (0i) cos (02)), we have 

(cos(0i)cos(02)) = 0.5(1 + V2) Jet = .000425, (19) 

and 

(sm(0i)sm(02)) = 0.5(1 - V2) Jet + CME = .000825. (20) 
The Jet = .0007692 and CME = -.00048 which is bigger than the RHIC value of -.00026. 

3.2 TMC at the LHC 

We need to add Transverse Momentum Conservation (TMC) to the LHC Pb-Pb same charge 
sign analysis. At first let us assume that scaling is not the same as for jets and could differ 
from V2. We can rewrite equations 19 and 20 as 

(cos(0i)cos(02)) = 0.5(1 + V2) Jet - 0.5(1 + ScaleT)TMC = .000425, (21) 

and 



(sm(0i)sm(02)) = 0.5(1 - V2)Jet - 0.5(1 - ScaleT)TMC + CME = .000825. (22) 
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Figure 4: The centrality dependence of the integrated three-particle correlator (cos(0q + 
0/3 — 2\E'rp)) (7) measured with four independent methods: the cumulant method using 
TPC tracks; the event plane estimation using the TPC, using the ZDC, or using the VZERO 
detectors. The correlation between same charge pairs and opposite charge pairs are indicted 
on the figure. 
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Figure 5: The centrality dependence of the integrated two-particle correlator (cos(0„ — (p^)) 
(6) measured by ALICE (close symbols) and STAR (open symbols). The correlation between 
same charge pairs and opposite charge pairs are indicted on the figure. 



The TMC value can be written in terms of the Jet value using equation 21. 

_ 0.5(1 + V2) Jet - .000425 
™^ - 0.5(1 + ScaleT) " ^^^^ 

TMC is a negative correlation and in equation 21 and 22 this negative sign has been 
taken into account making the TMC term of these equations positive. Prom equation 23 if 

TMC = 0.0 then Jet = .0007692. If we increase Jet then TMC will increase. How large is 
the Jet term in the LHC data? In the next subsection we use the opposite charge sign pairs 
in order to estimate the value of Jet and TMC. 



3.3 Opposite Sign Data and the TMC 

We can use the opposite charge sign data for centrality range 50% to 60% from Figure 4 and 
Figure 5. The value for (cos(0i + 02)) (7) is 0.0 and {cos{(j)i — 02)) (S) is .0033. We break up 
the 7 and 5 into an in plane term (cos(0i)cos(02)) and an out of plane term {sin{(t)i)sin{(f)2)) 
and obtain 

(cos(0i)cos(02)) = (1 + V2) Jet - 0.5(1 + V2)TMC = .00165. (24) 

Note that for ScaleT we have used V2 because for the opposite sign pairs there is no local 
suppression of near side pairs which increases the negative TMC effect between back to back 
particles. Also the 0.5 factor in front of the Jet term is gone. This doubling of the Jet term 
come from the fact that in p-p the opposite sign correlation for jets is twice the same sign 
correlation. 



3.4 LHC Predictions 

In the last two subsections we have determined two equations and two unknowns, namely 
the Jet term and the TMC term in equation 23 and equation 24. First let us consider the 
case where V2 is .105 and ScaleT = V2. From these two equations one gets Jet = .002217 
and TMC — .001448. When we substitute these values into equation 22 the CME value is 
-.00048. This value is the same as determined in the the first LHC subsection where there 
was no TMC term. 

What would an analysis of event shape give at the LHC? Like in the STAR event shape 
analysis we use half of the tracks of the TPC in order to determine the reaction plane. 
The other tracks are used to calculate (cos(0i + 4>2 — rp)) (7) for same sign. Also 
{cos2{(f) — \l'_Rp)) {V2) is also calculated for these same tracks. The V2 measures the shape of 
each subevent. For positive V2 we have more tracks in plane, while for negative V2 we will 
have more tracks out of plane. The value for Jet plus TMC in the correlator -7 should be 
zero when V2 is zero. The value of Jet plus TMC should be positive when V2 is positive and 
negative when V2 is negative. We have plotted Figure 3 again adding a solid line which is 
proportional to -7 that would be measured at ALICE if ScaleT is equal to V2 (Figure 6). 
The slope of the solid line is negative instead of positive like it is at RHIC. The difference in 
sign comes from the positive value of (cos(0i — (j)2)) (5) at LHC while at RHIC it is negative. 



The effect of tlie TMC term in tlie 7 correlator was determined at RHIC and was four 
times larger than V2. The same suppression of near side same sign charge pairs which is 
created by the strong chromo-electric fields inside the glasma flux tubes generated by the 
color glass initial state of the Pb-Pb collisions maybe at work. This suppression could greatly 
increases the negative TMC effect between back to back particles of the same charge sign. 
We do not know how strong this effect is at the LHC, but if we use the same value of ScaleT 
measured at RHIC ScaleT = .2836, then Jet = .002040 and TMC = .001094. When we 
substitute these values into equation 22 the CME value is -.000304. Again we can predict 
on Figure 6 (this time a dashed line) what should be measured at ALICE if ScaleT is equal 
to .2836. The slope of the dashed line has the same value as the solid line except it has the 
opposite sign. This slope is not as extreme as at RHIC. The slope is flatter than at RHIC 
because of a cancellation by the Jet term scaled down by ^2 (.105) with the TMC term scaled 
by ScaleT (.2836). At RHIC there is no Jet term. 



4 Summary and Discussion 

We have shown in this manuscript that The Chiral Magnetic Effect (CME) is consistent 
with the correlation analysis for same sign charge pairs done by the STAR experiment on 
mid-peripheral Au-Au Collisions yj/vjv=200 GeV. There are backgrounds present which can 
give signals that make the CME measurement hard to interpret. The STAR experiment has 
made a measurement that makes it possible to separate out the background from the signal 
for the same sign charge pairs. For the same sign charge pairs the Transverse Momentum 
Conservation (TMC) term is the only major background and it will very with the shape of 
particle production characterized by {cos2{(j) — ^ rp)) {v2). The V2 measures the shape of 
each event. For positive V2 we have more tracks in plane, while for negative V2 we will have 
more tracks out of plane. The value for TMC in the correlator (cos(0i + 02 ~ rp)) (t) 
should be zero when V2 is zero. The value should be negative when V2 is positive and positive 
when V2 is negative. Thus through this event shape analysis one can measure the CME. 

We have made a flt to same charge pairs correlations at a centrality value of 50-60% using 
equations 12-14. We have determined that at this centrality for Au-Au Collisions y^J/vjv=200 
GeV the value of CME is -.0002625 +.0000638 -.0000845. The TMC is two and half times 
bigger and of the same sign. The CME is a flow of charge out of the reaction driven by 
the magnetic fleld. This flow causes a charge separation and should be measurable. We can 
deflne in an event by event manner an in plane and out of plane charge separation 

AQiN = (N^; - N!i) - (N^^ - N^!) (25) 

and 

AQouT = (N^ - Nl) - (N^ - N^). (26) 

If there is a flow of charge out of plane then AQqut should be wider than the in plane AQin. 
For centrality 40-50% we show the AQ distributions (see Figure 7). In this Figure we also 
show that the rms difference between out of plane and in plane divided by the average rms 
of the two distributions. The fact that this difference is positive and not consistent with any 
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Figure 6: The STAR charge correlation difference (Aud) — {-^lr) and (A+A_)ud — A+A^)lr 
which is proportional to —7 for the same sign pairs (squares) and opposite sign pairs (dots) 
under the selection of soft particles low-pt {pt < 2 GeV/c) is ploted as a function of azimuthal i 
anisotropy shape of the sub-events. The solid line is a prediction of what would be measured 1 
at the LHC if the ScaleT is equal to V2 (.105). The dash is a prediction of what would be 
measured at the LHC if the ScaleT is equal to what it was at RHIC (.2836). ■ 
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Figure 7: Top panel: sample AQ distributions for 40-50% Au-Au a/s/vw = 200 GeV collisions 
in plane and out of plane. Bottom panel: is the rms difference between out of plane and in 
plane divided by the averge rms of the two distributions. 



known background other than CME, we conclude that the CME is the best and simplest 
explanation of the data. 

With the success of the RHIC Au-Au ^/sNN=200 GeV CME measurement we considered 
the LHC Pb-Pb y^SNN=2.76 TeV measurements. First we find that (cos(^i — ^2)) (<^) as 
measured by ALICE for the same charge sign pairs has changed sign. We assume that this 
is because of the importance of hard scattering or jets playing a greater role. By adding a 
jet term we are able to determine the amount of CME for the same 50-60% centrality. The 
value of CME is a little less than twice the value measured at RHIC -.00048 to -.00030. We 
can also make a prediction for the shape dependence of the event shape analysis if it is done 
at ALICE using their TPC (see Figure 6). This event shape analysis is a very important test 
of our approach to the same charge sign correlation data. However if the TMC effects scales 
like it did at RHIC the sign of the event shape correlation slope sign will be the same as at 



RHIC but of a smaller value. If this prediction is false then our analysis is false. One should 
note the lower range of the CME value falls with in the range of the RHIC result. One must 
also make the AQ measurements. The width of AQqut must continue to be larger than 
AQiN. 
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